The influence of clumping on predicted star wind parameters 
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We study the influence of clumping on the predicted wind structure of 0-type stars. For this 
purpose we artificially include clumping into our stationary wind models. When the clumps 
are assumed to be optically thin, the radiative line force increases compared to corresponding 
unclumped models, with a similar effect on either the mass-loss rate or the terminal velocity 
(depending on the onset of clumping). Optically thick clumps, alternatively, might be able to 
decrease the radiative force. 



1 Introduction 

Theoretical models of hot star winds led to the con- 
clusion that mass-loss has a significant impact on the 
stellar evolution in the upper HRD. This conclusion 
was supported by a relatively good agreement be- 
tween these models and observational data. 

However, both theoretical predictions and obser- 
vational values were derived, in a first approxima- 
tion, by assuming the stellar wind to be a smooth, 
spherically symmetric outflow. This picture might 
be inadequate due to the existence of a strong in- 
stability related to radiative line driving (see Feld- 
meier, this volume, and references therein). The 
influence of corresponding spatially inhomogeneous 
wind structures ( "clumps" ) on the basic wind prop- 
erties (i.e., mass- loss rate and velocity field) and the 
emergent spectrum was believed to be insignificant. 

This changed with the application of NLTE mod- 
els that were able to account for the influence of 
clumping on hot star wind spectra. The reason that 
the effect of clumping is not immediately apparent in 
the wind spectra is that most diagnostical features 
depend on the product yfClM^ where the "clumping 
factor" Cc > 1 relates the density inside the clump 



with the mean wind density (p) , 

P+ = Cc(p). 



(1) 



Consequently, spectra from winds with a large 
clumping factor but small mass-loss rate can mimic 
those from winds with weak clumping but large 
mass- loss rate. If Cc > 1, then the mass- loss rates 
derived from such diagnostics are overestimated by a 
factor of y/TTc. Fortunately, some spectral properties 
may be used to break this degeneracy, and to "obser- 
vationally" estimate the value of Cc. For example, 
Martins et al. (2005b) derived clumping factors of 
10 and 100 for studied Galactic 0-stars, implying a 
decrease of the estimated mass-loss rate by factors 



of 3 and 10 (see also Bourct, this volume, and Puis 
ct al., this volume, and references therein). 

Thus, if the actual mass-loss rates of hot stars are 
really lower than assumed before, then there is a sig- 
nificant discrepancy between the observations and 
theory. The reason for this discrepancy is not yet 
known. A possibility might be that clumping may 
shift the ionization balance in such a way that the 
final radiative force is lower, resulting in lower mass- 
loss rates. To test this possibility we have included 
clumping into our wind models and studied its in- 
fluence on the basic wind properties of hot stars. 



2 Model stars 

The wind models are calculated for 0-typc stellar 
parameters based on the recent calibrations by Mar- 
tins ct al. (2005a, see also Fig. 1). 



6.5 
6.0 



Schaller et al. (1992) 

85 Mq Martins et al. (2005) + 




45000 40000 35000 30000 25000 
7"ef[K] 

Figure 1: Parameters of studied stars in the HR di- 
agram. Ovcrplotted are the evolutionary 
tracks from Schaller et al. (1992). 
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3 Wind models 

For our study we used the spherically symmetric, 
stationary wind models developed by Krticka & 
Kubat (2004). These models solve the equations of 
statistical equilibrium together with the equations 
of radiative transfer. The calculated occupation 
numbers are used to derive the radiative force (in 
the Sobolev approximation) and the radiative heat- 
ing/cooling terms. This enables us to obtain the 
radial stratification of velocity, density and temper- 
ature in the wind and finally to predict the wind 
mass- loss rate M and terminal velocity v^o ■ For the 
studied 0-stars. M derived from our models depends 
on the stellar luminosity L and the effective temper- 
ature Toff on average as 




The advantage of our models is the self-consistent 
solution of the momentum equation, though the ra- 
diative transfer is treated in a simplified way. 

4 Optically thin clumps 

The assumption of optically thin clumps is widely 
used for studying the influence of clumping on the 
wind spectra (e.g., Martins et al. 2005b, Puis et al, 
this volume). 

To include optically thin clumps into our mod- 
els we modified the equations of statistical equilib- 
rium, by using an electron density p+ = Cc(po), 
opacity (x) = X^/C'c: and emissivity (rj) = ry+/Cc. 
The superscript -I- denotes values inside the (homo- 
geneous) clumps and the quantities inside brackets 
corresponding volume averages. 

4.1 Influence of clumping 

To investigate the influence of clumping on the stel- 
lar wind we have calculated a wind model of an O- 
type giant at Toff = 35 000 K, assuming the wind to 
be smooth (Cc — 1) close to the star (r < 2i?*), 
and to be clumped (Cc = 10) in the outer regions 
(r > 2i?H., see Figs. 2, 3). 

The presence of clumping leads to an increase 
of the electron density inside the clumps. Conse- 
quently, the recombination rates become higher and 
the wind less ionized (see Fig. 2). Since lower ions 
are able to accelerate the stellar wind more efficiently 
than the higher ones (due to a larger number of driv- 
ing lines), the radiative force increases, which, in 
our case, leads to an increase in wind velocity (see 
Fig. 3). 
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Figure 2: Influence of clumping on the ionization 
fractions of selected ions. 




Figure 3: Influence of clumping on the wind velocity. 



4.2 Radially constant clumping factor 

The influence of clumping on the wind parameters 
depends on the radial onset of clumping. If clump- 
ing starts above the critical point (below which the 
mass- loss rate is determined), then the terminal ve- 
locity increases (see Fig. 3). On the other hand, if 
clumping starts below the critical point, then the 
wind mass-loss rate becomes larger. 



Table 1: Average increase of M for constant Cc 



Cc 


1 3.16 


10 


31.6 


100 


Af(Cc)/A'/(Cc = 1) 


1 1.48 


2.15 


3.17 


4.57 



In case of radially constant clumping, the mass- 
loss rate increases significantly, and the predicted 
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wind-momcntuni rate is much higher than for a 
smooth wind at same parameters (see Fig. 4). 



5 Clumps larger than Lsob 

Individual clumps may be larger than the Sobolcv 
length I/Sob- Assuming the velocity gradient inside 
the clumps to be the same as in the corresponding 
smooth wind, we account for clumps being larger 
than Lsob by using the same modifications as in the 
case of optically thin clumps, but additionally in- 
creasing the Sobolev optical depth by Cc, and de- 
creasing the radiative force per unit of volume by a 
factor of Cr- 
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Figure 4: Influence of wind inhomogeneities on the 
modified wind-momentum rate. 



The inclusion of clumps larger than Lsob into our 
wind models leads to a decrease of the radiative 
force. Thus, and for radially constant Cc, the mass- 
loss rate decreases significantly (Fig. 4). 



Muijrcs (this volume): Indeed, wind inhomogeneities 
influence the predicted wind parameters. An ap- 
proach which roughly corresponds to the results of 
time-dependent simulations (i.e., clumps which are 
optically thin at most frequencies, "starting" above 
the critical point) does not improve the agreement 
between theory and observations. Clumps assumed 
to be larger than the Sobolev length (and start- 
ing below the critical point) may provide a better 
agreement between theory and observations, both in 
terms of mass-loss rates (Fig. 4) and P V ionization 
fractions (Fig. 5). Note, however, that the approach 
chosen by us strongly contrasts an important aspect 
of time-dependent simulations, namely that the ve- 
locity gradient inside the clumps is predicted to be 
much shallower than that of a corresponding smooth 
wind. This problem will be considered in a follow-up 
investigation. 
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Figure 5: Influence of different types of wind inho- 
mogeneities on the P V ionization. 
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6 Influence of porosity 

Wind porosity (Owocki et al. 2004) can be intro- 
duced into the wind models by additional decreasing 
the continuum opacity, Xeff = i+r^'(x) ' '^^ere h'r is 
the porosity length. The continuum cmissivity has 
to be modified by the same amount. 

This effect leads to a significant increase of the 
wind ionization. Thus, porosity leads to a decrease 
of the radiative force, and also the mass-loss rate 
may decrease if the wind is porous below the critical 
point (see Fig. 4 for the results for radially constant 
Cc and h'). 



7 Discussion 

Our results are in agreement with those derived from 
a similar investigation conducted by de Koter & 
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